What is the meaning of the statistical hadronization 
model? 



Francesco Becattini 

Universita di Firenze and INFN Sezione di Firenze, Via G. Sansone 1, F50019, Sesto F.no 
(Firenze), Italy 

E-mail: becattini@fi.infn.it 

Abstract. The statistical model of hadronization succeeds in reproducing particle abundances 
and transverse momentum spectra in high energy collisions of elementary particles as well as of 
heavy ions. Despite its apparent success, the interpretation of these results is controversial and 
the validity of the approach very often questioned. In this paper, we would like to summarize 
the whole issue by first outlining a basic formulation of the model and then comment on the 
main criticisms and different kinds of interpretations, with special emphasis on the so-called 
"phase space dominance" . While the ultimate answer to the question why the statistical model 
works should certainly be pursued, we stress that it is a priority to confirm or disprove the 
fundamental scheme of the statistical model by performing some detailed tests on the rates of 
exclusive channels at lower energy. 



1. Introduction 

The statistical model is a model of hadronization, thus aiming at reproducing the quantitative 
features of this process. Its founding ideas date back to Fermi and Hagedorn [2], though 
a basic and precise formulation of this model has been lacking ever since; providing such a 
formulation is among the goals of the present work. 

In the statistical hadronization model (SHM), the physical picture of a high energy collision 
is that of a QCD-driven dynamical process eventually giving rise to the formation of extended 
massive objects (called clusters or fireballs) which decay into hadrons in a purely statistical 
fashion. The number, as well as the kinematical and internal quantum properties of these objects 
are determined by the previous dynamical process and are thus not predictable within the SHM 
itself; they can be hopefully calculated with perturbative QCD, like in other hadronization 
cluster models [3] . A distinctive feature of the statistical model in comparison with other cluster 
models is that clusters have a finite spacial extension. This is actually a crucial assumption, the 
one which ultimately allows to make calculations. 

Probably, the best known model with relativistic extended massive objects is the bag model 
|I] and indeed the SHM can be considered as a model for the strong decays of bags. On the 
other hand, on the experimental side, there is now strong evidence of the finite extension of 
the hadron emitting sources in high energy collisions, from the observed quantum interference 
effects in the production of identical particles. It is therefore reasonable to take a finite volume 
of the hadron sources as a key ingredient for a hadronization model. 



In this paper, we first expound a precise formulation of the model starting from the very basic 
assumption of local statistical equilibrium (Sect. 2). We then comment on various criticisms and 
interpretations of the successes of this model, especially in elementary collisions (Sect. 3). The 
paper is concluded with a discussion about the possible fundamental physical meaning of the 
model (Sect. 4). 

2. The statistical model: a fundamentalist approach 

The basic idea of the model is very simple and it lies in two assumptions. The first is that 
in the late stage of a high energy collision, some extended massive objects, defined as clusters, 
are produced which decay into hadrons at a critical value of energy density or some other 
relevant parameter. The second, fundamental, assumption of the statistical model is that all 
multihadronic states within the cluster compatible with its quantum numbers are equally likely. 
This makes the model predictive as to the production rates of hadrons and resonances from the 
clusters so that SHM might be regarded as an effective model of the decays of hadronic extended 
relativistic massive objects. This idea was essentially introduced within the statistical bootstrap 
model by Hagedorn [2J, who, by identifying clusters with massive resonances, predicted the 
hadronic mass spectrum to rise exponentially. This seems to be still a very succesful prediction 
jS], but is neither implied nor required by the SHM alone: in principle, clusters need not to be 
identified with actual resonances to decay statistically. 

Despite the apparent simplicity of the key assumption of the SHM, it is not as straightforward 
as it might seem at first sight to calculate the cluster decay rates into different multi-hadronic 
states. These difficulties arise from the fact that the basic postulate only tells us that localized 
states are equiprobable, yet these states are essentially different from the observable asymptotic 
states. As we will show, such difference is not an issue when the volume is sufficiently large and 
can be disregarded in most applications where the canonical or grand-canonical ensemble are 
used, but it is relevant at a fundamental level of description and must be taken into account 
when the volume is small, i.e. less than O(10) fm 3 . 

Suppose that we can describe the cluster as a mixture of localized multi-hadronic states \hy) 
and, according to the basic assumption, all states have the same statistical weight. Then, one 
can write down a microcanonical partition function as: 



where Pj is the projector over all conserved quantities in strong interactions, namely energy- 
momentum, angular momentum, parity, isospin etc. It must be emphasized that these states \hy) 
are not the asymptotic observable free states of the Fock space which can be labelled with particle 
multiplicities for each species {N±, N2, ■ ■ ■ , Nk} = {Nj}, their momenta and helicities. Thus, 
the probability of observing a set of particles with four-momenta pi, ■ ■ ■ ,pn is not (hy\Pi\hy} 
and, moreover, it cannot be obtained unambiguously from the Eq. (|T|). 

In order to define a suitable probability of observing an asymptotic multi-hadronic state | /) , 
one can recast the microcanonical partition function Q by using the completeness of states 
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where Py = J2h v \hy)(hy\ is the projector on the localized states. We note that the last 
expression of in Eq. (J2J is a proper trace, whereas it was not in Eq. Q as the states \hy) do 



not form a complete set, i.e. they are not a basis of the Hilbert space. Looking at Eq. P]). it is 
tempting to set the probability pf of the state |/) as proportional to (/|PiPy|/)- Yet, one could 
have worked out f2 differently from Eq. for instance (if \hy) are properly normalized): 

n = J2( h v\ p i\hv) = J2( h v\ p vPi\h v ) 
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whence pf could now be set as proportional to (/IPyPjPyl/)- The latter expression is different 
from (/|PiPy|/) unlike [Pj, Py] = 0, which is not the case as long as Pj includes conservation of 
energy and momentum (we will see this more in detail in the next section). 
Which of the two probabilities: 

(f\P V PiPv\f) and (f\PiP V \f) (4) 

is the correct one? A well-defined probability should meet two requirements: 

• positivity; 

• respect conservation laws, i.e. pf = if |/) has not the same quantum numbers as the 
initial state; in other words, pf = if Pj|/} = 0. 

The leftmost expression in Eq. fulfills the positivity requirement in that: 

(/|PyPiPy|/) = (Pv/|Pi|Pv/> = a-^Py/IPflPy/) = o-^PiPv/IPiPy/) > (5) 

where we have used the hermiticity of Pj,Py and the fact that Pf = aPj through a positive 
divergent constant a . However, the second requirement is not fulfilled: pf turns out to be not 
vanishing even for states \f) which do not have the same energy-momentum as the initial state. 

On the other hand, the rightmost expression in (£Q) manifestly fulfills the conservation 
requirement because Pj has |/) as its argument, but, on the other hand, it is not positive 
definite. Positivity is recovered by changing the rightmost expression in Eq. (@J into: 

Pf <x {f\PiP v Pi\f) (6) 

i.e. by plugging one more projection operator on the initial state. Thereby: 

(/IPiPyPil/) = (Pi/|Pv|P*/) = (P*/|Py|Pi/) = (PvPif\PvPif) > (7) 

where now we have used the idempotency of Py. The definition © of the probability leads 
to a microcanonical partition function which differs from the proper one Q just by a positive 
(divergent) constant which is anyhow irrelevant for the calculation of averages. Indeed: 

tr(PiPyPi) = tr(P?Py) = atr(PiPy) = afi (8) 

where Eq. © has been used. 

The leftmost definition of the probability in Eq. @ was in fact used in ref. [S] to work out 
the rates of multi-hadronic exclusive channels {Nj} in the SHM. Although breaking energy- 
momentum conservation, the expressions of these rates are the same as those obtained from 

The divergence of this constant is owing to the non-compactness of the Poincare group. This can be understood 
by considering the projector on energy-momentum S 4 (P — P op ) 



integrating Eq. © over momenta of final particles. An attractive feature of this definition is 
that it can be written formally, with Pj = as: 

p/<x|</|Pv|;>i 2 (9) 

i.e. it looks similar to a transition probability. However, the expression © seems to be the best 
suited because it naturally meets both aforementioned basic requirements. 

According to the definition (JBJ), the cluster is regarded as the mixture of states: 

Y,Pi\hv)(h V \Pi (10) 

hy 

unlike in the leftmost definition in Eq. (|1J). where the mixture turns out to be: 

Y,\bv){hv\Pi\h v )(hv\ (11) 

hv 

We think that a mixed state where Pj|/iy) are equiprobable, like in (|lUj) is the most appropriate 
definition of microcanonical ensemble because any state in the mixture has actually the same 
quantum numbers, including energy and momentum, of the cluster. 



2.1. The microcanonical ensemble 

Once a satisfactory definition of the probability of observing an asymptotic multi-hadronic state 
is found, that is Eq. @, we can start working it out to obtain tractable formulae. We start by 
first developing the projector Pj defining the microcanonical ensemble. 

In principle, in the microcanonical ensemble, all conserved quantities should be included: 
energy, momentum, angular momentum, parity, internal charges and C-parity (if the cluster is 
neutral). The correct way to implement these conservation laws is to project the multi-particle 
states onto the irreducible state of the full symmetry group which defines the initial state, i.e. 
the hadronizing cluster. The full symmetry group is the product of the extendend Poincare 
group 10(1, 3)1, the isospin SU(2), the U(l)'s related to conserved additive charges and the 
discrete group Z2 of charge conjugation if the initial state is neutral. Accordingly, the projector 
Pj can be factorized as: 

Pi = Pp,J,A,- P X P /,/aPQ (12) 

where P is the four-momentum of the cluster, J its spin, A its helicity, tt its parity, \ its C- 
parity, I and 1% its isospin and its third component and Q = (Qx, ■ ■ ■ ,Qm) a set of M abelian 
(i.e. additive) charges such as baryon number, strangeness, electric charge etc. Of course, the 
projection P x makes sense only if I = and Q = 0; in this case, P x commutes with all other 
projectors. 

The projector Ppjx% onto the irreducible state (transforming according to an irreducible 
unitary representation v of 10(1,3)1) with definite four-momentum, spin, helicity and parity can 
be written by using the normalized invariant measure \x of the Poncare group as: 

Pp,J,x,, = ~ E dimi/ / M9z) D»\g z )\ U(g z ) (13) 
z z=\,n J 

where z is the identity or space inversion g z G 10(1, 3)^., D l '{g z ) is the matrix of the irreducible 
representation v the initial state i belongs to, and U(g z ) is the unitary representation of g z in 
the Hilbert space. Similar integral expressions can be written for the projectors onto internal 
charges, for the groups SU(2) (isospin) and U(l) (for additive charges). Although projection 



operators cannot be rigorously defined for non-compact groups, such as Poincare group, we 
will maintain this naming relaxing mathematical rigour. In fact, for non compact-groups, the 
projection operators cannot be properly normalized so as to P 2 = P and this is indeed related 
to the fact that \i) has infinite norm. Still, we will not be concerned with such drawbacks 
thereafter, whilst it will be favourable to keep the projector formalism. Working in the rest 
frame of the cluster, with P = (M, 0), the matrix element D v \g z )\ vanishes unless the Lorentz 
transformations are pure rotations and this implies the reduction of the integration in ()13|) from 
IO(l,3) T to the sub group T(4) (g> SU(2) <8> 0- I n fact, the general transformation of the 
extended Poincare group g z may be factorized as: 

g z = T(x)ZA = T(x)ZL ft (0R (14) 

where T(x) is a translation by the four- vector x, Z = l,n is either the identity or the space 
inversion and A = !_„(£) R is a general orthocronous Lorentz transformation written as the 
product of a boost of hyperbolic angle ^ along the space-like axis n and a rotation R depending 
on three Euler angles. Thus Eq (jlMjl becomes: 

Pp,J,x,* = \j2^^Jd i xjdAD^J(x)ZA)rU(J(x)ZA) 

= ^ E ^ / d ^/ dAe lP W(A)» U(J(x))U(Z)U(A) (15) 

where z = if Z = I and z = 1 if Z = In the above equation, the invariant measure d 4 x 
of the translation subgroup has been normalized with a coefficient l/(27r) 4 in order to yield a 
Dirac delta, as shown later. Furthermore, dA is meant to be the invariant normalized measure 
of the Lorentz group, which can be written as [7]: 

dA = dl_ n (OdR = sinh 2 £d£^dR (16) 

dR being the well known invariant measure of SU(2) group, £ G [0, +oo) and f2„ are the angular 
coordinates of the vector n. 

If the initial state \i) has vanishing momentum, i.e. P = (M, 0), then the Lorentz 
transformation A must not involve any non-trivial boost transformation with £ ^ for the 
matrix element D v (l\) 1 * not to vanish. Therefore A reduces to the rotation R and we can write: 

Pp,j,x,, = \ E 7^4 / d ' x ( 2 J + 1) / dR e lP ' W J (R)>- U(J(x)) U(Z) U(R) (17) 



z=i 



Since [Z, R] = 0, we can move the U(Z) operator to the right of U(R) and recast the above 
equation as: 

Pp,J,x,n = jty!^* elP ' XU ( J ( x ))( 2J + 1) / dR ^ J ( R )a* U(R) ' + 7r ^ (n) 

= 5\P - P op )(2J + 1) J dR D J (R) X X * U(R) l+^iHl (18) 



The Eq. (|18|) is indeed the final general expression of the projector defining the proper 
microcanonical ensemble with P = (M, 0), in which all conservation laws related to space- 
time symmetries are taken into account. The appeal of the above expressionresides in the 



factorization of projection operators onto the energy-momentum P, spin-helicity J, A and parity 
7r of the cluster. 

Also the projectors onto isospin Pij 3 and onto additive charges Pq in Eq. (fT2]l can be given 
an integral expression by using the invariant SU(2) and U(l) group measures. The projector 
on a state with definite C-parity x can be simply written as (I + xQ/2 where C is the charge- 
conjugation operator. 

In most calculations, conservation of angular momentum, isospin, parity and C-parity is 
disregarded and only energy-momentum and abelian charges conservation is enforced. 2 This is 
expected to be an appropriate approximation in high energy collisions, where many clusters are 
formed and the neglected constraints should not play a significant role. On the other hand, they 
are important in very small hadronizing systems (e.g. pp at rest (Hj) whereby the full projection 
operation in Eq. Q12|) should implemented. For the restricted microcanonical ensemble, it can 
be easily seen from Eqs. (|12I18|) that the projector can be rewritten as: 

Pi = S\P-P op )S Q)Qop (19) 

and, if |/) is an eigenstate of four-momentum and charges, Eq. as: 

Pf <xu5\P-P f )5 Q)Qf (f\P v \f) (20) 

where w is a divergent constant, i.e. the whole space-time volume. 

2.2. The effects of finite volume 

The second part of the calculation involves the projection onto localized states 
Py can be written as: 

{Nj},k 

where {Nj} are the occupation numbers of the particles in the cluster and k the variables 
labelling their kinematical modes, e.g. three integers in case of a parallelepipedon box with 
fixed or periodic boundary conditions. It should be pointed here that in the usual statistical 
model calculations, the interactions are taken into account by including all known resonances 
as free particles (in the cluster) with a distributed mass, according to a formalism developed 
by Dashen, Ma and Bernstein 0^3- I n principle, the use of H21|) to calculate probabilities 
like in © entails some difficulty because a localized state \{Nj},k) is not an eigenstate of the 
actual particle number, which is defined in terms of the operators creating and destroying free 
asymptotic states over the whole space. In fact, a A-pion state in the cluster has non vanishing 
components on all free states of the pion field, i.e. on the states with 0, 1, 2, . . . pions. Therefore, 
the projection should be performed in a full quantum relativistic field approach by identifying 
localized states as states of the quantum fields associated to particles and vanishing out of the 
cluster region. Hence, the projector Py should be rather written as, in case of only one scalar 
particle: 

Py = / D^)^\ (22) 
Jv 

and Eq. © developed accordingly. In Eq. (|2*2*)l \tp) = ® x |V'(x)) and Dip is the functional measure; 
the functional integration must be performed over all functions having as support the cluster 
region V . Altogether, determining production rates involves the calculation of the statistical 
mechanics of a field in the microcanonical ensemble. 



The projector 
(21) 



2 Note that the relevant set of states is still defined microcanonical ensemble and we will comply with this 
tradition. 



Nevertheless, if the cluster size is sufficiently larger than the Compton wavelenght of the 
particles involved, quantum field corrections are expected to be small and the eigenstates of 
particle number operators in the whole space essentially correspond to those of the particle 
number operators in the cluster. Otherwise stated, even though a localized iV-pion state has 
non vanishing components on all free states of the pion field, the dominant one will be that on 
the asymptotic N pion states. The largest implied Compton wavelenght in a multi-hadronic 
system is indeed the pion's one A^ ~ 1.4 fm; this is the minimal size of the cluster below which 
quantum field corrections cannot be disregarded. 

Hence, for clusters which are sufficiently larger than A^ we can use the approximation: 

({N 3 },p\{N 3 },k) + iff {Nj} = {Nj} (23) 

where p labels the set of kimatical variables (namely momenta and helicities) for the asymptotic 
free states |/) = \{Nj},p). Now, by using Eq. (|2*H|) the probability of a final state lf2Ti|) can be 
calculated and reads: 

p {Njhp oc 5\P - »)<Vv «,. EKW'flW'^l 2 ( 24 ) 

i k 

where the p^s are the four-momenta of the particles in the final state. The rightmost factor 
in the above equation can be calculated as a cluster decomposition and, in the framework of 
non-relativistic quantum mechanics, a relevant expression has been obtained in ref. in the 
limit of large volumes and in ref. |S] taking into account the finite volume. If j labels the hadron 
species j = 1, . . . , K and p now the set of particles' four- momenta: 

p W >, - «v - £») n e {T1)N3 N +Hj( l! i+1)Hj 3 ( 25 ) 

i 3 {hnj} Un-=l n j ', 1 

where {h n } is a partition of the integer Nj in the multiplicity representation , i.e. Nj = 
En/=l n iS ; H i = Dn/=l h n 3 and: 

^-UjWJv***^ ^ (26) 

are integrals over the cluster region V, q being the cyclic permutation of the integers 1, . . . , n\. 
For large volumes, the dominant term in the cluster decomposition (|25[) is obtained by taking 
{hfiA = (Nj,0, . . . , 0), implying c; = I and reads: 

i 

3. Phase space dominance, Lagrange multipliers and all that 

Despite its apparent success in reproducing observables related to hadronization process like 
particle multiplicities and transverse momentum spectra j!21 1141 115j the statistical model 
is not very popular among high energy physicists. Besides the fact that, thus far, the featured 
apparent statistical equilibrium is not derivable from QCD (and it will remain so for a probably 
long time), one of the most bothering point seems to be the presence of thermodynamical 
quantities and chiefly temperature. Moreover, this fairly constant hadronization temperature 
(i.e. around 160 MeV |15| ) is amazingly close to the estimated critical temperature of QCD and 
this obviously raises the question of its meaning. 
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V N i{2Jj + l 
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In the basic micro canonical formulation of the model, described in the previous sections, 
one deals with mass and volume of clusters, but it can be easily shown that if those become 
sufficiently large one can perform calculations in the canonical ensemble, which is far easier to 
handle, thereby introducing temperature through a saddle-point expansion [^JEI]- ^ n the- case of 
hadron gas, this is possible at relatively low values of masses and volumes ^QlEl, around 8 GeV 
and 20 fm 3 . Yet, talking about temperature in such small systems seems to be daring for the 
received wisdom of most physicists who, as soon as the word temperature is spoken, are led to 
think of a large system which has undergone a long cooking process before reaching equilibrium. 
It is widely believed (the author is included) that this cannot occur after hadronization at a 
level of formed hadrons through inelastic collisions: the system expands too quickly to allow 
this. If statistical equilibrium is genuine, it must be an inherent property of hadronization itself, 
i.e. hadrons are born at equilibrium as stated by Hagedorn many years ago jlUj and reaffirmed 
by others more recently [T31 ITH1 ITTH 1271] . 

Therefore, there have been some attempts to account for the success of the statistical model 
whose conclusions may be roughly clustered as follows: 

(i) the results of the statistical model can be obtained from other models with some 
supplementary assumption or invoking some special, so far neglected, mechanism; 

(ii) the statistical model grasps some truth of the hadronization process, but the apparent 
thermal-like features are an effect of a special property of the quantum dynamics governing 
hadronization, which tends to evenly populates all final states: this is defined as phase space 
dominance; 

(iii) the results of the statistical model are somehow trivial, due to the large multiplicities 
involved which eventually make the multi-hadronic phase space almost evenly populated. 

In the following I will comment on specific papers discussing this subject, whose attitude, 
on the basis of my personal understanding, is assigned to one (or more) of the previous points. 
I apologize in advance with the quoted authors for possible misunderstanding and too limited 
summary of their thought. 

Ideas of the class (i) are proposed e.g. in refs. The typical exponential shape of 

the thermal spectra are explained in the framework of the string model, by adding to the basic 
picture additional fluctuations of the string tension parameter k. The effect of the fluctuations 
is to broaden the gaussian shape of the pt spectra in the string model, turning it into an 
exponential one. Of course, this mechanism is one of the possible choices of nature, though 
very difficult to disprove. In general, it is certainly possible to make an existing model more 
complicated to account for some otherwise more straightforward result in another model, and 
this is precisely where the criterium called Occam razor intervenes: between two models equally 
able to explain observations, the most economical should prevail. It is fair to say here that 
the string model has been tested against more observables and that the SHM should be tested 
against the same set of observables. Still, it is also true that the effective implementations of the 
string models are plagued by the need of many free parameters to reproduce the data and this 
raises many doubts about its predictive power [23]. In fact, there is an ongoing work I24j 
to implement SHM as hadronization model in an event generator to allow testing observables 
other than multiplicities and single particle inclusive transverse momentum spectra. 

The paper by Hormuzdiar et al [2E] is the one where the idea (ii) is certainly argued more in 
detail. The basis of the whole argument is the similarity between the (classical) phase space of 
the set of particles {Nj}, obtained by integrating (|2T|) over momenta: 

j 4 (p*-E^) ( 28 ) 

i 




where N = YljNj, and the general expression of the decay rate into the channel {Nj} of a 
massive particle (cluster) in relativistic quantum mechanics: 
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where | Mu | 2 is the Lorentz-invariant dynamical matrix element governing the decay. Assuming, 
for sake of simplicity, all spinless particles, |M/-j| 2 may in principle depend on all relativistic 
invariants formed out of the four-momenta of the N particles, as well as on all possible isoscalars 
formed out of the isovector operators I,. Suppose that |Mjj| 2 = a , so that the whole dynamics 
reduces to introduce the same multiplicative constant a for each particle in the channel. Then, 
it is possible to calculate quite easily the generating function of the multi-particle multiplicity 
distribution starting from Eq. (|29j) : 



G(Xi,...,Xk) = E r {w,}II^. 
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(30) 



where a Fourier decomposition of the four-dimensional delta has been used and the integral Wick 
rotated by using z = \x. If P 2 is sufficiently large, we can expand the above integral around the 
saddle-point zq obtained by solving the equation: 



^ + 5z^(2vr)3 / 9,. 6 
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If P = (M, 0) it is not difficult to realize that zq = (/?, 0) and the generating function can be 
approximated as: 



G(Ai, . . . , X K ) ~ exp 



E ( 27r )3 / 9 ,. 
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so that the mean number of particles of the species j reads: 
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which is very similar to a thermal distribution: 
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(32) 



(33) 



(34) 



were not for the different measure in the momentum integral. As it should be clear from its 
derivation, the constant (3 in Eq. (|33|) is certainly not a temperature, rather a soft scale parameter 



which is related to the effective finite interaction range. Yet, the ratios of average multiplicities 
of particles of different species mimic a thermodynamic behaviour. This is the so-called phase 
space dominance. The authors of ref. [25] work out a more specific example based on QED and 
they conclude, quite reasonably, that a fairly good fit to particle multiplicities may be provided 
if integral expressions like (|33|) are used instead of an actual Boltzmann integral. I want to even 
reinforce their statement by adding that the actual fits to particle multiplicities in e + e~, pp 
and other collisions relied on supplementary assumptions which are not expected to be exact 
in a basic statistical model framework, so that deviations from the "pure" statistical model 
predictions may arise which can be of the same order of the difference between the actual SHM 
and the formula (|3*3l . 

Thus, I subscribe to the argument in ref. [13] but it should be stressed that the just described 
phase space dominance is a highly non-trivial assumption. In fact, the recovery of a thermal- 
like expression like ()33|) ought to a very special form of the matrix element |Mfj| 2 , where both 
the dependence on kinematical and isospin invariants was disregarded. If a different form, still 
perfectly legitimate and possible, is assumed, the thermal-like behaviour is spoiled. For instance, 
one could have: 



with a generic factor f(arrij)g(Ij) for each particle depending on its mass rrij and its isospin Ij 
and on a single scale a whose dimension is the inverse of an energy; the factors o?M in Eq. (|35j) 
are introduced in order to make the average particle multiplicities in the large multiplicity limit 
proportional to the mass of the cluster M: 



It is not difficult to realize that the production function Eq. (|36|) might be dramatically different 
from the thermal one. It should be emphasized that also a factorizable dynamical matrix element 
depending only on masses and isospins like in Eq. 1)35(1 is quite an exceptional one. In fact, in 
principle, there could be dependence on other independent invariants like (pi +Pj) 2 , Ij • lj, 
Ij • (Ij x Ifc) etc. Therefore, an observed phase space dominance in multihadron production is 
not a trivial fact and tells us something important about the characteristics of the underlying 
non perturbative QCD dynamics, besides providing us with an empirically good model. 

Similar arguments are presented in ref. [21)] and, more extensively, in ref. [2T] where the 
concept of phase space dominance is even more explicitely defined. There, quantities like the 
previous (5 arising from a saddle-point asymptotic expansion and mimicking a temperature are 
called "Lagrange multipliers" just to emphasize the difference from an actual temperature. 
However, in the effort of analyzing the meaning of the statistical model results, two very 
questionable statements are introduced: 

• that the so-called Lagrange multipliers have no physical meaning even for a properly defined 
phase space integral like (|28|): 

• that the phase space dominance is trivial when the average multiplicities are very large. 

For the second point, the counter-argument is straightforward: just take f(m) = j4 exp(Cm 2 ) 
(though odd it might look) or, alternatively, g(I) = AI + C in Eq. (|3l)|) with A, C positive 
constants depending on centre-of-mass energy and the thermal shape of mass production function 
is destroyed for any multiplicity. 

The first point is more subtle and requires a somewhat general discussion because there 
seems to be some confusion as to what deserves to be called "thermal" and, conversely, what is 
only "statistical" . If the word "statistical" is used to to mean some property of the dynamical 




(35) 




(36) 



matrix element of being independent of most kinematical variables, like that leading to Eq. (j33[) . 
then of course it has nothing to do with a proper thermal thing. If, on the other hand, the 
word "statistical" means, like in SHM, equal probability in phase space, where phase space is 
appropriately measured with d 3 xd 3 p for any particle like in Eq. (|28[) . and a volume is involved, 
then "statistical" and "thermal" can be taken as synonimous (for purists only for sufficiently 
large volumes) because there is no quantitative difference between them. In fact, what makes 
the difference between Eq. (|33|) and a proper thermal formula is the measure in momentum 
space and the absence of a volume. If, in a proper statistical mechanical framework, the two 
conditions of statistical equilibrium and finite volume are met, temperature can be defined (e.g. 
through a saddle point expansion) no matter how the system got to statistical equilibrium and 
even in absence of an external bath. Many authors (e.g. [28 ) take the definition T" 1 = dS/dE 
where S is the entropy, a well defined quantity for any closed system. All other definitions of 
temperature, be a Lagrange multiplier for the maximization of entropy at a fixed energy |29j . 
or a saddle point of the microcanonical partition function, should converge to the same value in 
the limit of large volumes and are therefore physically meaningful temperature. Macroscopically 
inspired definitions requiring physical exchange of energy with a heat reservoir are too restrictive, 
and certainly not suitable for heavy ion collisions as well, where such a heat reservoir does not 
exist. On the other hand, these definitions must coincide with the most general definition based 
on statistical mechanics. 

In the same spirit, some authors |29| try to make clear a distinction between the temperature 
determined in the SHM by fitting particle abundances and a "proper" temperature which would 
be achieved through inelastic reinteractions of formed particles. The former is called Lagrange 
multiplier for the maximization of entropy, just to emphasize the difference. Again, I would 
like to stress that there is no actual quantitative difference between those two temperatures 
so that a hadronization temperature, if confirmed, can be properly called a temperature. One 
can certainly make a distinction as to how statistical equilibrium was achieved, which is as 
important as the statistical equilibrium itself, but if energy is equally shared among all possible 
states within a finite (possibly large) volume, temperature is temperature no matter how the 
system got to statistical equilibrium. What would make the exponential parameter fitted in the 
framework of SHM different from an actual temperature can be only a quantitative difference, 
like e.g. the difference between (3 in Eq. (|33j) and (3 in Eq. (|3^|) . 

What can be done then to distinguish between a genuine statistical-thermal model and 
other possible pseudo-statistical models like the one leading to the formula 1)33(1 ? Besides 
kinematical features, it would be desirable to bring out effects related to the finite volume, which 
is a peculiarity of the statistical model. Indeed, the study of average inclusive multiplicities 
or inclusive pt spectra does not allow clearcut conclusions because those observables are 
not sensitive enough to different integration measures (i.e. Vd 5 p versus d 3 p/2e) and much 
information is integrated away. A much more effective test would be studying the rates of 
exclusive channels, i.e. F{ N .}/T{ N /y, which are much more sensitive to the integration measure 
in the momentum integrals and the shape of dynamical matrix element. Unfortunately, exclusive 
channels can be measured only at low energy (some GeV) where none of the conservation laws, 
including angular momentum, parity and isospin, can be neglected, as pointed out in ref. [S] 
where pp annihilation at rest has been studied in this framework. This makes calculations rather 
cumbersome and difficult from the numerical point of view. None of the numerous previous 
studies in literature has tackled the problem without introducing approximations unavoidably 
implying large errors in the calculations. Fully microcanonical calculations including both four- 
momentum and angular momentum conservation have not ever been done, and only recently 
the increased computing power and purposely designed techniques allowed the calculation of 
averages in the microcanonical ensemble, yet only with energy and momentum conservation 

HE3EEE]. 



4. What is the meaning of it? 

Now that we have discussed in some detail the foundations of the statistical model, and possible 
interpretations of its success, we are finally left with an inevitable question: what is the meaning 
of this model in the framework of the basic theory of strong interactions, QCD? Otherwise stated, 
is it possible to show from a more fundamental theory that extended massive objects such as 
clusters exist and that the statistical filling of their multihadronic phase space effectively occurs? 
Or, alternatively, that QCD implies a similar phenomenon (though quantitatively distinct and 
distinguishable), called phase space dominance? As yet, we are not able to answer this question 
because QCD has not been solved in the non-perturbative regime. Therefore, we will try to 
argue about some simpler issue. 

A first issue is the meaning of the mixture of states (jlUj) that we have used to describe cluster 
decays. From a quantum mechanical viewpoint, a mixture of states is only a mean to describe 
our ignorance of the state of the system, which is always supposed to be a pure one, be it 
entangled or not. We do not want here to slip into fundamental quantum mechanics problems 
like decoherence and measurement, which may render a mixture of states an objective description 
of the system. Just to make this issue a concrete one in our perspective, it suffices to mention a 
(low energy) collision creating one cluster: of course this should be described with a pure state. 

Let \i) the pure quantum state of a cluster; we can instance think of this state as that which 
can be calculated in the bag model in terms of free parton fields states confined within a finite 
region. We can write the transition amplitude to a localized multi-hadronic state within the 
cluster: 

(h v \T\i) oc (h v \TPi\i) = {P t hv\T\i) (37) 

where the last equality follows from the conservation laws, that is the transition operator T 
depends on the hamiltonian of strong interactions and ought to commute with the projector 
onto conserved quantities. We can build up a basis of the Hilbert space including the \hy) 
vectors by adding to them the multihadronic states localized outside V, i.e. the region denoted 
with V. We can then write: 

i = ]T|M<M+EM<M (38) 

hy ky 

Essentially, the results of the statistical model can be recovered by assuming: 

(hy\T\i)=0 y\hy) 

\(Pih v \T\i)\ 2 = \c hv \ 2 = C (39) 

where C is a constant, independent of the state Pj|/iy). The first of the two equations in (|39j) 
states that no transition can occur to a state outside the cluster volume; the second, that the 
transition probability is uniform for all localized states with the same quantum numbers as 
the cluster itself. In a sense, these statements amount to restate the Hagedorn's hypothesis of 
a resonance as being made of a uniform superposition of hadrons and resonances. From the 
previous assumptions and using Eqs. (|37j ) .(|38 |) . one can calculate the transition amplitude to an 
asymptotic state |/): 

</TO = (/l (EIM<V| + EIVXV|) T\i) = J2(f\hv)(Pih V \T\i)=J2(f\ h v)ch v (40) 

\hy hy ) hy hy 

so that: 

\(f\T\i)\ 2 = \J2(mhv}c hv \ 2 = J2\(f\^ h v)\ 2 C + E (f\Pihv){Pih! v \f)c hv c* h , v (41) 

hy hy hyj^h'y 



The first term in the right hand side of above equation is just proportional to @. So, the 
statistical model results are fully recovered if: 



or, in other words, if the amplitudes Ch v = \[C exp(i0/ lv ) defined in 1)39(1 have random phases 
(f>h v , so to make the cross-term sum vanishing. 

Hence, we have actually rephrased the question whether the statistical model can be an 
effective model for the hadronization process actually driven by QCD on the question whether 
the conditions (|39|) and (|42jl apply in a QCD-inspired picture. Since it is presently not possible 
to answer to this question either, we are left with the more approachable problem of verifying 
the predictions of the statistical model more thouroughly, as we have discussed at the end of 
previous section. 

5. Conclusions 

We have discussed in some detail the ideas and the interpretations of the success of the statistical 
model in reproducing soft observables in high energy collisions. It is certainly crucial to 
understand the why of this success from firts QCD principles, but in the meantime it is useful 
to stick to a more pragmatic attitude and ask ourselves whether we can test this model more 
deeply than what has been done as yet. Particularly, by testing the model against exclusive 
channel rates, we can assess whether the thermal- like features of inclusive particle production 
show up at high energy because of the quasi-independence of dynamical matrix elements in 
the soft non-perturbative regime (phase space dominance). In fact, it is difficult to bring out 
deviations from a genuine statistical model from the analysis of inclusive quantities only because 
too much information is integrated away. On the other hand, such deviations should show 
up in more detailed observables, like, e.g., exclusive channel rates. In this regard, relevant 
data are available only at low energy (some GeV in centre-of-mass frame) and this requires 
the implementation of full micro canonical calculations, which have never been done without 
introducing too drastic approximations. We have outlined an appropriate framework (in Sect. 2) 
for the full micro canonical formulation of the model, on the basis of group projection techniques. 
This the first step to implement the calculation; numerical work is currently ongoing. 
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